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Abstract 

On the basis of the full analytical solution of the overall unitary dynamics, the 
time evolution of entanglement is studied in a simple bipartite model system evolving 
unitarily from a pure initial state. The system consists of two particles in one spacial 
dimension bound by harmonic forces and having its free center of mass initially lo- 
calized in space in a minimum uncertainty wave packet. The existence of such initial 
states in which the bound particles are not entangled is pointed out. The entanglement 
of the two particles is shown to be independent of the wavepacket mean momentum, 
and to increase monotonically in a time scale distinct from that of the spreading of 
the center of mass wavepacket. 

1 Introduction 

Since the early development of Quantum Mechanics, a number of questions have been 
raised concerning its interpretation and status facing the physical description of reality. 
'Entanglement' is one such question, pointed out and developed by Schrodinger pQ in 
the wake of the paper by Einstein, Podolski and Rosen on the completeness of Quantum 
Mechanics [2]. The situation considered by Schrodinger in that context hinges on the 
description of the quantum mechanical state of each one of two separate parts of a com- 
posite system, after they have interacted for some time. He observes that even if, prior 
to the interaction, the state of each of the constituent parts is described by its own quan- 
tum state vector (or, in the words of Schrodinger, by its p si- function), in which case the 
quantum state of the composite system is described by their product, this is in general no 
longer possible after interaction has run its course, the constituent parts having become 
'entangled' by the intervening interaction process. This entanglement is thus associated to 
the non factorizability of the quantum state of the composite system after the interaction, 
a trait of quantum mechanics that "enforces its entire departure from classical lines of 
thought" pQ. 

After the breakthrough of Bell [3], which opened the way to experimental probing [1] 
into quantum peculiarities relating to entanglement [5], this concept has become a 
subject of intense research [6], particularly in connection with the fact that it provides 
crucial resources for the new research domains of quantum information and quantum 
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computation [8]. Prominent issues in this context have been, on the one side, devel- 
oping the characterization and quantification of entanglement, especially as applicable to 
quantum states which are themselves not represented in terms of a state vector, or a 
psi- function [6l [9]; and, on the other side, understanding the dynamics of entanglement 
specially in realistically open systems, i.e. systems undergoing interactions with a com- 
plex 'environment' [9]. Here one aims particularly at gaining control over the degrading 
effects of external interactions on the valuable entanglement properties residing in the 
system of interest. These effects are generally referred to as 'decoherence' processes, and 
have been given also a positive role in connection with accounting for the practical elu- 
siveness of in principle allowed quantum superpositions of macroscopically distinguishable 
states |10| 111], i.e., in connection with the 'construction of a classical looking world' on 
the basis of quantum mechanics. 

More recently, investigations of the effect of decoherence processes on the time evolu- 
tion of entanglement present in a quantum system led to the finding of the rather surprising 
effect which became known as ESD, "entanglement sudden death" (or, more euphemisti- 
cally, "early stage disentanglement") [12] . which consists in the complete demise, within 
a finite time, of the entanglement initially present in the system of interest. This effect 
has in fact been observed in terms of a phenomenological effective, non unitary quantum 
dynamics which was simulated with photon polarizations in the role of two level sys- 
tems [13]. Subsequently ESD has been also found to be present in extremely simple closed 
(tripartite) systems evolving unitarily (in fact, multiperiodically) from initially entangled 
psi- functions |14j . showing that environment induced decoherence is not a requisite for 
this particular phenomenon in the dynamics of entanglement. Furthermore, the fact that 
one of the subsystems involved in the ESD process does not interact with the remain- 
ing parts relativises the role of dynamical correlation processes in the time evolution of 
entanglement. 

This work deals with the dynamics of entanglement in the even simpler context of 
a bipartite model system which is completely soluble in analytical terms, using tools 
available in any introductory course on quantum mechanics. The system consists of two 
particles in one spacial dimension bound by harmonic forces and having its free center 
of mass initially localized in space in a wave packet of minimum uncertainty. By using 
center of mass and relative motion dynamical variables, the complete unitary quantum 
dynamical evolution of this state can be easily written down explicitly. Assuming that 
the two particles are in a stationary state (e.g., the ground state, for simplicity) of the 
relative motion, the non stationary character of the initial state reduces therefore just to 
the dispersive time evolution of the center of mass wavepacket. While the center of mass 
and relative 'subsystems' are initially not entangled, having their own psi-functions, and 
remain so at all subsequent times, the same is in general not true when one considers each 
of the two individual particles as subsystems. It turns out, in fact, that the two particles 
are in general entangled, and that their entanglement is time dependent as a result of 
the dispersive time evolution of the center of mass wavepacket. While the presence of 
entanglement is of no surprise, in view of the fact that the two particles are strongly 
correlated dynamically by being in a bound state of the harmonic interparticle force, 
one finds nevertheless that the initial center of mass wavepacket can be prepared so that 
the two particles are actually not entangled initially. The dispersive time evolution of 
the center of mass wavepacket is found to induce a monotonic increase in time of the 
entanglement of the two particles, albeit in a time scale distinct from that associated with 
the dispersive spreading of the center of mass wavepacket. These facts exemplify in very 
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simple terms the surprises and subtleties one may find in the relation of entanglement of 
subsystems and its dynamics to the unitary dynamical evolution of the system as a whole. 

The paper is organized as follows. In section 2 the complete characterization of the 
model system is given and its dynamics is fully worked out analytically, both in terms of 
center of mass and relative coordinates and in terms of particle coordinates. The pres- 
ence of entanglement in the latter description is also discussed. Section 3 constitutes the 
main body of the work. It is divided in three sub-sections, dealing respectively with the 
description of the quantum state of entangled subsystems, with some pertinent quanti- 
tative measures of entanglement and with a rather detailed discussion of the dynamics 
of entanglement in the model system. Section 4 contains some concluding remarks. The 
explicitation of a more technical point relating to the adopted measure of entanglement is 
given in an Appendix, for completeness. 



2 Alternate descriptions of a flying two-body bound state 

We consider a system consisting of two (spinless) particles in just one dimension, which in- 
teract through a translationally invariant attractive potential. If the masses and positions 
of the particles are {mi,Xi} and {771,2, £2} respectively and V{x 2 — x\) is the potential 
describing the interaction of the two particles, the dynamics of the system is described by 
the Hamiltonian 

where p\ and p 2 are the momenta canonically conjugated to x\ and x 2 - Note that the 
interaction potential V correlates the dynamics of the two particles. A common simplifying 
practice in this context is to deal with this correlated two-body problem in terms of an 
alternate set of dynamical variables describing the center of mass and relative motions 
rather than in terms of the particle dynamical variables themselves. The reason for this 
is that in terms of the new variables one has only to deal with two independent, albeit 
effective, one-body problems. In fact, the appropriate transformations to center of mass 
variables X, P and relative variables x,p are 

R = ■ , P = Pl +p 2 ; r = x 2 -xi, p= ■ , (2.1) 

777,1 + fni mi + "72 

in terms of which the the Hamiltonian of the two-body problem appears in the guise 

Here the new mass parameters are, respectively, the total mass M and the reduced mass 
fi of the particles, which in terms of the particle masses can be expressed as 

A/f . mi777 2 



mi + 777 2 



This reduction procedure can be used both in classical and in quantum contexts. In the 
former case it involves a canonical transformation, while the latter it involves a unitary 
transformation of the basic canonical variables. Even though this is in fact rather trivially 
evident, it is worth noting explicitly that in both cases the system maintains a composite 
character, but which in the latter rendering reduces to two non-interacting parts. 
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We now recall the simple quantum mechanical treatment of the system, taking as a 
specific instance of the two-body interaction potential a simple harmonic potential V(x) = 
/MjJ 2 x 2 /2. The stationary solutions of the Schrodinger equation for the composite system 
are well known and simple to obtain. They consist in the product of a center of mass 
momentum eigenfunction with a stationary state of the harmonic oscillator equation which 
governs the dynamics of the relative motion. Taking the relevant state to be the oscillator 
ground state one has 

TT^Vb V ^ 

Note that one faces here the well known normalization problem in the case of the con- 
tinuous spectrum of the center of mass Hamiltonian. We will not have to deal explicitly 
with this, however, since we consider a (non-stationary) normalized gaussian wavepacket 
of width B and mean momentum HK as the wavefunction describing the initial state of 
the center of mass part. The relative part remains in its chosen ground stationary state, 
so that our normalized non-stationary initial state can be written as 

$(#,r) = — e -^e lKM — . (2.2) 

In order to obtain the time evolution of such a globally non-stationary state it is most 
convenient to have it rewritten in the momentum representation by means of a double 
Fourier transform 

= J e-^ R e-^ r ^(R,r)dRdr = 4^ e ~2B^-Kf e -&? . ( 2 . 3 ) 

Since the relative state is an oscillator eigenstate, the Fourier transform involving the 
relative coordinate is strictly not necessary for this purpose, but it will be convenient 
further on for reverting the description to the original particle variables, instead of the 
center of mass and relative variables. 

Time evolution can now be implemented by just bringing in the appropriate time- 
dependent phase factors associated with the different center of mass momentum eigen- 
functions as well as with the relative oscillator eigenfunction. This leads to the well 
known translation and spreading of the center of mass wave packet, with no other effect 
than an irrelevant overall time-dependent phase factor originating from the internal energy 
eigenstate. More formally, the time evolution is represented by the action of the evolution 
operator 

U(t,0) = <g> e ~^ e1 *, 

which reduces to numerical phase factors when acting on the momentum space represen- 
tation. Thus, after elapsed time t, the momentum space wave-function (|2.3p becomes 

V 7T 

The simplicity of this result is basically due to the independence of the center of 
mass and relative dynamics. In particular, the stationary character of the relative mo- 
tion state implies that its contribution to the joint time-evolution reduces to an overall 
time dependent phase factor. The time evolution of the center of mass wave packet is 
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itself more involved, as different momentum components acquire different time-dependent 
phases which govern its dispersive displacement. The center of mass and relative motion 
parts of the composite system have however their independence fully uncompromised by 
these time-dependences, as clearly exhibited by the factorization of the overall wavefunc- 
tion into center of mass and relative motion factors at all times. 

Now let us return to our former decomposition of the system, in terms of its two 
constituent particles. In momentum representation, this can be done on the basis of the 
solution written in eq. (|2.4p . through reverting to the dynamical variables relating to the 
individual particles. The appropriate transformation here is (cf. eqs. (|2.ip ) 

K - jk + fc, £ - mifc2 ~ m2fcl 

where k\ and k 2 are the wave numbers associated with the momentum eigenfunctions of 
particles 1 and 2 respectively. In terms of these, the time evolved state appears as 

— (ki+k 2 ) 2 p -2B 2 (k 1 +k 2 -K) 2 p -if* -f^(mifc>-m2fcl) 2 




(ifel, k 2 ,t) = A\ — e W l*l+**J g-^ fcl + fc2 -Kr e - lt t e -^mi« 2 -m 2Kl , _ (2 _ 5) 



An important feature of this alternate rendering of the state of the two-particle system 
is that the dependence on k\ and on k 2 does not factor for all times t > 0, as can be seen 
by the occurrence of terms involving the product k\k 2 in the exponents. Note however 
that factorization will in fact occur at t = for a particular choice of width parameter 
B of the center of mass wavepacket. In fact, when this parameter is tuned to the value 
satisfying 

B* Bl = = ^ ,2,, 

these non- factoring terms drop out, and one obtains a product wavefunction also in terms 
of the individual particle variables. This is clearly a result particular to our choice of a real 
gaussian envelope function for the center of mass wavepacket, combined with the gaussian 
form of the relative ground state in the harmonic particle-particle interaction potential. 

The factorizability of the wavefunction fj2.4|) at all times t allows for associating the 
quantum state of the center of mass and relative parts of the system with a respective 
factor of that wavefunction. Each of these two parts is therefore endowed with its own 
independent quantum mechanical probability amplitude. Non factorizability when it is 
expressed in terms of particle variables as in eq. (|2.5|) . on the other hand, implies that 
the joint probability amplitude includes such correlation involving the two particles which 
impedes an independent description of each one of them in terms of probability amplitudes. 
This is just what has been called by Schrodinger (in 1935, [I]) the quantum 'entanglement' 
of the two particles. 

Our findings so far indicate moreover that a definite distinction should be made between 
entanglement and correlations between parts of the composite system. In fact, in the state 
described by eqs. (|2.4p and (j2.5fl . the two particles are, at all times, and in particular at 
t = 0, strongly correlated by finding themselves bound in the stationary ground state of 
relative motion. The initial state factorization in terms of particle variables for B = Bq 
provides thus for a counter example disallowing direct association of such correlation with 
entanglement. The factorized amplitudes of the wavefunction cast in terms of the particle 
variables, eq. (|2.5p . in a generic factorized instance B 2 = fib 2 /M at t = 0, appear as 

ipi(ki) = Nie M m i 1 M and ^2(^2) = N 2 e M 2 M , 
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Ni and N 2 being the appropriate normalization constants. It is apparent in these expres- 
sions that, while the quantum kinematical property of amplitude factorization ensures lack 
of entanglement, each one of the two factor amplitudes involve parameters referring to the 
two particles and bearing on their dynamical correlation. 



3 Quantum states of entangled particles 

According to the rules of quantum mechanics, knowledge of the joint wavefunction (|2.5p 
provides us with the most complete information one can obtain about the system itself. 
It must then, in particular, also provide us with the most complete information one can 
obtain about each of its constituent particles, considered separately. The formal tool 
usually employed to extract these more restricted pieces of information from the amplitude 
(|2.5p describing the quantum state of the two-particle system involves first setting up the 
so called density matrix, defined as 

pikuh-k^k'^t) = ^(k 1 ,k 2 ,t)^*(k' 1 ,k 2 ,t). (3.7) 

The 'diagonal' part of this, i.e., the part with k\ = k[ and k 2 = k' 2 , corresponds just to the 
joint probability distribution (in momentum space) of the two particle system, at time t. 
Moreover, from the normalization of the wavefunction ^?(k\, k 2 , t) it follows that 

J dk'i J dk'ip{k x M\KX^)p{KX-AA^) =p(ki,k 2 ;k' 1 ,k' 2 ;t). 

The double integration works here just as a matrix multiplication, so that this relation in 
fact shows that, considered as representing an operator, the density matrix (|3.7p associ- 
ated with the normalized wavefunction (|2.5p is hermitian and idempotent, and therefore a 
projection operator. 

Now the mean value of any observable pertaining to the first particle can be 
obtained in terms of the density matrix as 



{O w ) t = I dhj dk[0^ ki J dk 2 p{k u k 2 ; k[,k 2 ; t) 
dh [ dk[0$ ' pW(fci,fci;t), 



where is the momentum representation of the observable under consideration, 

Since this works for any observable pertaining to particle 1, it follows that the object 



p W (h, k[;t) = / dk 2 p(ki,k 2 ;k[,k 2 ;t), 



known as the reduced density matrix for particle 1 contains all information concerning 
this particle, when considered independently of the rest of the system. The operation 
leading to the reduced density matrix is usually referred to as taking a partial trace of 
the density matrix p(k\, k 2 ; k[, k 2 \ t) over particle 2, and amounts to an all inclusive ac- 
ceptance of whatever properties it would have correlated to the particle under scrutiny. 
A similar argument allows one to collect all information concerning particle 2 which is 
similarly inclusive with respect to properties of particle 1, by taking a partial trace of 
p(ki, k 2 ; k[, k 2 ;t) over this particle, i.e. 
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p {2) {k 2 ,k' 2 ]t) = J dki p(k 1 ,k 2 ;k 1 ,k 2 ;t). 
Note that, again in view of the normalization of the wavefunction (|2,5p . one has 

J dk 1 pW(k 1 ,kr,t) = J dk 2 pW{k 2 ,k 2 ;t) = 1. (3.8) 

An application of these ideas to the specific case of the two-particle system as consid- 
ered here goes as follows. First, the density matrix reads explicitly (cf. eqs. (|2,5p and 



p{k u k' x M,k' 2 ,t) = — e ^(^+^) 2 -«+ fe i) 2 )) e - 2B2 (^+^-^) 2 +( fc i+ fc i-^) 2 ) x 

7T 

xe -j^{(rnik2-m2k 1 ) 2 +(m 1 k' 2 -rn 2 k' 1 ) 2 ) /g 



Next, the reduced density matrix for particle 1 can be obtained by taking a partial trace 
of this over particle 2, which amounts here to setting k 2 = k' 2 and performing a gaussian 
integration over the momentum variable k 2 . This gives the somewhat lengthy result 



^(k 1 ,k[;t)= (3.10) 



8BbM e ~ M2s2+m i 62 ., , (fci-fc') 2 , j . r2 , 1 , m (M f ~ 9 MB !'o (fci+fcH 2 

o 16(M 2 B 2 +m 2 b 2 ) y 11 ^4K B £ (k\-k\) g M 2 B 2 +m 2 b 2 y X 

TTy/M 2 B 2 + m'(b 2 



x e 



2KB^ M f B -^ 2 ) (fa+fci) , (M 2 B2+m , 62)(fc 2 +fcf) -if ( M2 ^. fc . ) (fc?-fcf ) 



In view of the complete symmetry of the density matrix fj3.9[) under an interchange of 
indices lf>2, the reduced density matrix for particle 2 can be immediately read off from 
eq. ()3.10p simply by replacing all indices 1 by 2 and vice versa. 



3.1 Reduced density factorization and entanglement 

An important feature of the reduced density matrix ()3.10j) and also of its companion 
p( 2 \k 2 ,k' 2 ;t) is that their dependence on primed and unprimed variables does not factor 
in general. This is in fact clear in view of the exponentials involving {ki ± k^) 2 . Non- 
factorizability implies that the reduced densities cannot be obtained from single quantum 
amplitudes in the same way that the full density matrix (|3.7p is obtained from the wave- 
function (|2.5p . and indicates the entanglement of the two particle system in the selected 
quantum stat^E The reduced densities are however hermitian objects in the sense that 

p^*(k i ,k' i ;t)=p^(k' l ,k i ;t), i = l,2 . 

A useful representation of the reduced densities can therefore be obtained in terms of the 
solution of the eigenvalue problem 

1 Recall that factorization may occur in special cases, such as the one pointed out earlier (see eq. H2.6|l 'l in 
which the wavefunction (|2.5|l factors, namely t = 0, and B — > Bo = b^J p,/M. In this particular situation it 
is possible to associate probability amplitudes to each of the two particles, although each of the amplitudes 
will involve the mass parameters of both particles, in addition to 'shared' parameters such as the oscillator 
parameter b and the width Bo of the center of mass wavepacket. 



7 



Hermiticity ensures that the eigenvalues are real and the eigenfunctions can be made 
orthonormal. In addition, from the fact that the mother density (|3.7p is a projection oper- 
ator, one can actually deduce that the eigenvalues are in fact non negative, i.e. \n\t) > 0. 
The reduced density p^(ki, k[; t) can therefore be written in the form 

P ®M;t) =E A n ) W xttHk^x^iklt) (3.11) 

n 

which engages all the (orthonormal) eigenfunctions with nonvanishing eigenvalue. Note 
also that, in view of the normalization property (|3.8p . one has ^ A„ (t) = 1. 

A rather non trivial outcome of this representation is that the two reduced densities 
originating from a given mother density matrix of the form (|3.7p have the same eigenvalues, 
their eigenfunctions being related as 

X ^(k 2 ,t) = N n (t) J dhx^ih, t)^(h,k 2 ,t) 

where N n (t) is an appropriate normalization factor. Furthermore, the joint wavefunction 
itself may be written in terms of the reduced densities eigenfunctions as (see Appendix)) 

= Y,V^Xn\k 1 ,t)x£ ) (k 2l t), (3.12) 

n 

in what is usually referred to as the Schmidt decomposition [8]. 

Therefore, what one can do in general is to write the reduced densities as linear com- 
binations of orthogonal projection operators with non negative coefficients of unit sum, 
also known as convex sums of orthogonal projection operators. Note that cases in which 
dependence on primed and unprimed variables factors in the reduced densities are duly 
covered as cases in which a single nonvanishing, hence unit valued, eigenvalue occurs in 
eq. (|3.1ip . This makes the reduced densities themselves projection operators. 

Whenever the convex sum (|3.1ip involves more than a single term the same happens 
with eq. (|3.12j) . which reveals a connection of non factorization of the dependence on 
primed and unprimed variables in the reduced densities and the non factorization of the 
dependence on k\ and k 2 in the wavefunction of the composite system. This is precisely the 
technical characterization of entanglement given by Schrodinger in his paper of 1935 [1]. 
Thus, entanglement is revealed in this case by the fact that the reduced densities have more 
than a single nonvanishing eigenvalue (hence needing more than a single eigenfunction for 
its representation) . In view of eq. (|3.12p one can see this as scars left from the process of 
'orphaning' of quantum amplitude correlations implied by multiple terms in eq. (|3.12p . 

Finally, it is worth stressing that the preceding analysis has been performed on objects 
ultimately related to the wavefunction of the composite system at a given time t. Therefore 
all the derived quantities are in general time dependent in a way which can be ultimately 
traced back to the unitary time evolution of the wavefunction of the composite system. 
In particular, if one devises a quantitative measure of entanglement, one will be able to 
follow its evolution in time. 
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3.2 A simple quantitative measure of entanglement 



The coefficients in the convex sum over orthogonal projection operators in eq. (|3.1ip . given 
the condition ^ n X n (t) = 1 which follows from normalization of the wavefunction of the 
composite system, can be seen as a probability distribution over the set of projections onto 
the eigenvectors of the reduced density. Furthermore, as the occurrence of more than a 
single nonvanishing weight signals entanglement of the two parts of the considered system 
which are under consideration, one is led to the use of entropy-like functions as they are 
used in statistical mechanics to provide for quantification of the degree of entanglement. 
A common choice for that is the so called von Neumann entropy |15|. [6] , given as 

S[p (1) (t)] = -^A n logA n . 

n 

This is a function endowed both with additivity and convexity properties [16], but its 
evaluation is often more costly than alternate functions, such as the so called '2-Renyi' 
entropy S2R and the 'linear' entropy A [6], which are defined respectively as 



Sfcfl[p (1) (*)] = - log and A[p( 1 )(t)] = ^A„(l-A„) = l-^A,2. 

n n n 

The evaluation of the sum appearing in these objects can be carried out in a simple way 
in terms the reduced density matrix itself, by noting that 

]Ta2 =Tr[p( 1 )(t)] 2 = f dh f dk 2 p^(ki,k 2 ;t)pW(k 2 ,k 1 ;t). 

n J J 

An identical result is obtained if the reduced density p <y2 \k\,k2;t) is used instead of 
p W(k 1 ,k 2 ;t). 

Both of these functions are non-negative and have the desirable convexity properties, 
although the linear entropy is not additive. This fact is however of little relevance in this 
context, and since, unlike the values of S2R, the values of A are bounded above by 1, we 
will use the linear entropy for quantifying the entanglement of the two particles. 

Evaluation of the linear entropy for the reduced density obtained in eq. (|3.10p gives 



A(t) = l-Tr[pW(t)] 2 = 1 



tft 2 (M 2 B 2 +m 2 b 2 )(M 2 B 2 +m 2 b 2 ) 
16M' 2 B 2 b 2 "I M 4 B 2 b' 2 



b tr 1 , 

In the last expression, tr = MB 2 /his the characteristic time associated with the dispersive 
spreading of the center of mass wavepacket, and r has been introduced as the characteristic 
time associated with the time dependence of the linear entropy. It is given by 



TBI 



m 2 b 2 \ / mVa 1 \ 

1 + wm){ 1 + Mhp)- (3 - 14 > 

an expression which is symmetric under the interchange of indices 1 2, consistently with 
the fact that the same result for the linear entropy is obtained if one uses p^ 2 \k2,k 2 ,t), 
instead of p^ l \k\, k[, t), in its evaluation. 
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Figure 1: Overall view of the initial linear entropy landscape as a function of log 10 mi/m2 and B/b. 
The logarithmic scale for the mass ratio emphasizes the symmetry under interchange of m\ and 
m,2. Most interesting features are seen to develop along the plane mi/rri2 = 1 (log 10 m\/mi = 0), 
mainly around B/b — 0.5. 



3.3 Initial entanglement and its dynamics 

A first, general feature to be noted in connection with the expression obtained for the 
linear entropy A(f) is that it turns out to be independent of the mean center of mass 
momentum K. This indicates Galilean invariance of the dynamics of entanglement in the 
present case, and is similar to what happens in connection with the dispersive spreading 
of wavepackets in position space. Actually this feature can be examined in very simple 
terms by considering the action of the center of mass momentum translation operator 

G(C) = e< x = e< m 

on the representation of the wavefunction of the composite system given in eq. (|3.12j) . 
The fact that the two position operators x\ and x<i commute allows for the factorization 
of G(C) as Gi(C) ® G 2 (C), so that 

G({)*(ki,k2,t) =52VK(G 1 (OxP(ki,tj) (G 2 (C)xl 2) (A; 2 ,t)). 

n 

The center of mass momentum translation therefore does not affect the amplitudes \f\n 
(also known as the coefficients of the Schmidt decomposition) which govern the entan- 
glement. The eingenvalues of the reduced densities are therefore preserved, and so is, in 
particular, the value of the linear entropy. 

Consider next the value of the reduced entropy for the class of initial states of the form 
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Figure 2: Closer view of the most interesting section of the initial linear entropy landscape as 
a function of log 10 mi/m2 and B/b. The section shown here contains the point mx/m^ = 1, 
B/b = 0.5 where the valley of A(t = 0) at mi/mz = 1 for B/b > 0.5 branches into two valleys 
where A(0) = 0, separated by a ridge still at mi/ma = 1, for B/b < 0.5. 



(|2.3p . It is given in terms of the the mass parameters mi, m% and of the gaussian width 
parameters b, B as 



A = A(i = 0) = 1 



_6 tb 

S T 



_6 



1 



b 2 



1 + 



1 



m\b 2 
APB 2 



1 + 



B 2 (1 + mi/nfi) 2 



+ 1 



.(3.15) 



As shown, it can be expressed in terms of just the two dimensionless quantities B/b and 
mi/m,2- Their value, supplemented by one additional scale setting quantity, which can be 
conveniently taken to be the total mass M, completely determines particular initial states 
of the assumed form. Note that the linear entropy is independent of the scale setting, so 
that its behavior over different domains in the manifold of states of the form (|2.3p may 
be mapped by analyzing Ao as a function of B/b and m\/m2- An overall view of the 
behavior of this function can be seen as a surface plot in fig. [TJ 

We now discuss some features of this initial entanglement landscape which can be 
scrutinized by analyzing eq. (|3,15p m some detail. First of all, one easily finds that, for 
any given value of the mass ratio m\/m2, the relevant extremum of Ao as a function of 
B/b occurs for B/b —> ^m\mil (l+wii/n^). Substitution back into the expression for Ao 
shows that one has at these points Ao = 0. We thus re-obtain, in the guise of the vanishing 
of the initial value of the linear entropy, the parametric conditions for the factorization of 
the initial wavefunction when cast in terms of particle variables. 

All the possible values of B/b associated with this set of unentangled initial states lie 
in the interval < B/b < 0.5, as follows from the range available for mi/m2, namely 
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A(0), mi/m 2 =1 



0.02 




B/b 



0.425 0.45 0.475 



0.525 0.55 



Figure 3: A cut of the surface shown in fig. [5] at the plane a = m\/m,2 = 1, showing values of the 
maxima (for ft < 0.5) and minima (for ft > 0.5) of the initial linear entropy. The plotted range of 
/3 = B/b is the same as in fig. [2] 

< m\/m2 < oo. They are therefore confined to a front vertical slice of fig. [TJ In order 
to explore the behavior of the initial value of the linear entropy for B/b > 0.5 one may 
fix the value of this variable and study the dependence of the initial linear entropy on 
m\/m2- What one finds in this way is that there is an extremum at m\/m2 = 1 for all 
possible values of B/b, in addition to two other extrema with m\jm2 values given by 



which are real, non-negative numbers only for B/b < 0.5. It turns out that the extremum 
at mi/m,2 = 1 is in fact a (nonvanishing) minimum of the initial linear entropy for B/b > 
0.5, which however becomes a maximum for B/b < 0.5. In this range, the additional 
extrema are the minima for which A(0) = 0, which have already been identified earlier. 
The fact that there are two such minima is due to the symmetry of the linear entropy 
under interchange of the masses mi and 77t2- As a check on this, one can easily verify 
that (7711/7712)— (7711/7712)4. = 1 for all values of B/b. At B/b = 0.5, A(0) has a very flat 
minimum (a zero of the fourth order) as a function of 7711/777-2 at 7711/777,2 = 1. These 
features are displayed in fig. [2j which covers the relevant domain of the initial linear 
entropy landscape, namely 0.031 < 7711/777,2 < 31 and 0.42 < B/b < 0.55. A cut of this 
surface graph, showing the values of A(0) at mi/772,2 = 1 as a function of B/b is shown 
in fig. [3l In this graph the ranges B/b < 0.5 and B/b < 0.5 correspond respectively to a 
ridge and to a valley in the A(0) surface. 

The time dependence of the linear entropy measure of particle-particle entanglement 
consists in a simple, monotonic approach to the upper bound A = 1 as seen in eq. (|3.13p . 
It involves an algebraic expression which is of the same form as that which governs the 
dispersive spreading of the center of mass wavepacket albeit involving a time scale r ^ tb- 
Since tb is the relevant time scale in the unitary dynamical evolution of the initial state of 
the composite system, it is both natural and convenient to use this time scale to analyze the 
dependence of the characteristic time r of the dynamics of entanglement on the different 
initial states of the form ()2.3|) . 

The quantity t/tb given in eq. (|3.15p is seen to be closely related to the initial 
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Figure 4: Overall view of the t/tb landscape as a function of log 10 mi/m2 and B/b. Note the 
relation to the features of the corresponding landscape for the initial linear entropy. 

state value of the linear entropy. In particular, it is also independent of the scale setting 
parameter M, and can be written explicitly as 



This implies that the main general features of the A(0) landscape are essentially carried 
over to the t/tb landscape, as can in fact be clearly in fig. |H One finds here the same 
overall structure of valleys and ridges seen in the corresponding initial linear entropy plot, 
in fig. [TJ A similar comparison can be made for the close-ups into the more intricate 
region of the m\/mi X B/b plane, figs. [5]and[2j As a consequence of inversion and of an 
additional b/B factor, however, one sees that t/tb is bounded by unity below, while A(0) 
is similarly bounded above. For sufficiently small values of B/b the characteristic time for 
the evolution of entanglement becomes substantially longer than tb- 

Note however that more tightly localized wavepackets spread faster (have smaller values 
of tb), so that this does not imply very long 'absolute' times for entanglement evolution. 
Using the value of M for scale setting, the reduced mass \x and the value of b are set by the 
values of mi/m^ and of the oscillator frequency u, the wavepacket width B being finally 
fixed by B/b. The limiting value of the 'absolute' characteristic time r for tightly confined 
center of mass wavepackets fixed in this way is 




(3.16) 
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Figure 5: Closer view of the valley bifurcation at mi/ma = 1, B/b = 0.5 in the t/tb landscape. 
Compare with fig. [2J 

4 Concluding remarks 

We have examined in some detail the quantum entanglement of two particles in one dimen- 
sion whose center of mass evolves as a free gaussian wavepacket with mean momentum hK, 
when they are bound in the ground state of a two-body harmonic interaction potential. 
While the fact that the occurrence of entanglement is of no surprise in view of the strong 
mutual correlation (e.g. in position space) of the two particles [17], one can easily see in 
this case that the entanglement may nevertheless be made to vanish in the initial state by 
an appropriate choice of the width parameter of the center of mass wavepacket. The fact 
that particle-particle correlations can, at least in some cases, be represented in unenta- 
gled form (which, in the present case, means being represented by factorized amplitudes) 
indicates that the two concepts are in fact quite distinct. 

Given the stationary character of the bound state of the two particles, the dynamics 
of their entanglement is driven by the unitary, dispersive spreading of the center of mass 
wavepacket. In particular, it also comes out as being independent of the mean value of the 
center of mass momentum, a fact which can be understood in terms of the preservation 
of the coefficients in the Schmidt decomposition of the two particle wavefunction under 
center of mass boosts. The time dependence induced by the unitary driving dynamics 
on the particle-particle entanglement is algebraically akin to that governing the change 
in time of the width parameter of the center of mass wavepacket, but involves its own, 
different time scale. The two time scales become even more distinct as the initial width 
parameter of the center of mass wavepacket becomes appreciably smaller than the bound 
state size parameter. In fact, we have seen that the quadratic decrease, as B — > 0, of the 
characteristic time tb for wavepacket spreading is not followed by the characteristic time 
for the evolution of entanglement, which approaches a finite limit involving properties of 
the two particle bound state. 

One must note also that the monotonic increase in time of both the spreading of the 
initial width B of the center of mass wavepacket and of the linear entropy A used to 
measure the degree of particle-particle entanglement hinge on the particular choice made 
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for the initial state. A still particular, but different choice leading to also different behavior 
may be devised by taking advantage of the time-reversible character of the overall unitary 
dynamics. In fact, one may use as an alternate initial state ^>(k\, k 2 ; 0), the time reversed 
counterpart of the state resulting from our adopted initial state (|2.3p after it evolves for a 
time T, i.e. 

^(h,k 2 ;Q) = ^f*(-h,-k 2 ;T). 

Due to time-reversal invariance of the overall dynamics, the time evolution of this state 
is given as \P(A;i, k 2 ;t) = Vl/*(— ki, —k 2 ; T — t), so that the time evolution from ^f(ki, k 2 ; 0) 
to ^(ki, kT) is traced backwards in time. As a result, the center of mass wavepacket will 
shrink and the linear entropy will decrease as t increases from to T. In particular, if 
^(ki, k 2 ; 0) is a particle-factorizable initial state, the linear entropy decreases smoothly 
from a non zero value to zero in the finite time T, as in cases of the so called entanglement 
sudden death [12] , but immediately and also smoothly rebounds to the monotonic increase 
for subsequent times. 

Finally, it is worth stressing that our simple model system provides for instances of 
quantum states consisting of parts which are strongly correlated dynamically but which 
are nevertheless not entangled. This is in a way complementary to the situation treated 
by Schrddinger [1] in the wake of the Einstein, Podolski and Rosen [2] argument, namely 
one in which entanglement persists after interaction between the two parts has ceased. 
Both cases point to a picture in which entanglement is ultimately a quantum kinematical 
feature, even if circumstantially affected by quantum dynamical processes. 
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Appendix 

We give here, for completeness, a brief, elementary calculation showing that the two 
reduced densities pWQci, k 2 ), i = 1, 2 have the same eigenvalues and that the wavefunction 
ty(k\,k 2 ) can in fact be written in the form shown in eq. f|3.12j> . cf. ref. |15j . Chap. VI, 
section 2. We omit the time variable, which is not relevant here. 

Note first that the set of eigenfunctions of p^ can be extended to a complete or- 
thonormal set {%n (k±)} in the Hilbert space of particle 1 by including eigenfunctions 
with eigenvalue zero. Consider, in addition, a complete orthonormal set {u u (k 2 )} in the 
Hilbert space of particle 2. Then, for each of the Xn^ define 

X ( nHh) = ^U u (k 2 ){x { n ) U v ^), 
v 

the bracket being the scalar product in product space. This may be checked to be an 
eigenfunction of p^ 2 \ In fact, write this reduced density as 

P (2 \k 2 ,k' 2 ) = J2 S dk i S dk^^h^ikuh^ik^k'^HK) 

so that 



15 



n' v 

= A«xi 2) (fc 2 ). 

The sum over v in the last expression involves the reduced density p^ and is in fact 
just \n 8 n n'i and the integral over k% yields the function Xn\k 2 ). This shows then that 
this function is an eigenfunction of the reduced density p^ with the same eigenvalue as 
its parent eigenfunction Xn\k\) of pW. From now on we may thus omit the eigenvalue 
superscripts. 

(2) 

Next check the norms of the functions xh, ■ This is easily done by evaluating 

V 

t<2) (2) 

showing that the functions Xn (^2) = Xn (fej/vAn constitute an orthonormal set. On 
the other hand, k 2 ) may be expanded in the product base {Xn\k\)u u {k2)}'- 

^(h,k 2 ) = ^^xl 1) (^iK(fc2)(xl 1) ^,^) = 

n v 

= YXnHh)X { nHk2)=Yv^x£Hh)X { *\k2) 
n n 

which establishes eq. (|3.12p . 
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